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We suggest that box shape trap enables to observe intrinsic properties of the repulsive Hubbard 
model in a fixed doping in contrast to the conventional harmonic trap bringing about spatial vari- 
ations of atom density profiles. In order to predict atomic density profile under the box trap, we 
apply the directly-extended density-matrix renormalization group method to 4-leg repulsive Hub- 
bard model with the open boundary condition. Consequently, we find that stripe formation is 
universal in a low hole doping range and the stripe sensitively changes its structure with variations 
oiU/t and the doping rate. A remarkable change is that a stripe formed by a hole pair turns to one 
by a bi-hole pair when entering a limited strong U/t range. Furthermore, a systematic calculation 
reveals that the Hubbard model shows a change from the stripe to the Friedel like oscillation with 
increasing the doping rate. 

PACS numbers: 71.10.Fd, 71.10.Pm, 74.20.Mn, 03.75.Ss 



Ultra-cold atomic Fermi-gas has attracted not 

only atomic gas community but also several physicists 
studying strongly-correlated electron systems. The rea- 
son is that the so-called "optical lattice" [3] formed by 
counter laser beams creates a periodical lattice potential 
described by the tight-binding model and the "Feshbach 
resonance" enables to access to the Hubbard model with 
the repulsive on-site interaction. Thus, a research goal 
in the optical lattice with the Feshbach tuning is to di- 
rectly observe several controversial issues due to strong 
correlation as seen in High-Tc cuprate superconductors 
and other metal oxides in controllable manners [4]. 

Generally, ultra-cold Bose and Fermi gases are trapped 
inside a harmonic trap to avoid the free expansion of 
atoms. The optical lattice is created inside the trap po- 
tential by utilizing the field modulation in the standing 
waves of lasers. Then, the model Hamiltonian on Fermi 
atoms (e.g., in 1-dimensional (1-D) case) is described b; 
the Hubbard model with a harmonic trap potential [H, 
given by 

i^Hubbard = -t^{0.^ja^ic^ + ^.C.) + U ^nj^Tlii 



V 



N - 



I — 1 , (1) 



where, a two-component Fermi gas is assumed, a-^ is 
the creation operator of a Fermi atom with pseudo-spin 
(J =t or I, nicr{= a\^aia) is the density operator of the 
i-th site, N is the total number of sites, and the sum- 
mation ii^j) in the first term describing the tunneling 
between lattice site is usually taken over the nearest 
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neighbor sites. The on-site repulsive interaction U is 
controllable via Feshbach resonance, and the presence of 
the harmonic trap potential is characterized by the last 
term including V . In the model Hamiltonian ([1]), it is 
well-known that the atomic density profile gives spatial 
variation, whose typical one is composed of the central 
Mott domain n^^r ~ 1) and the periphery metal- 
lic edges i^^Tiicr < 1) in a large repulsive U/t range. 
Such a variation is universal irrespective of the space 
dimension as long as the harmonic well type of trap is 
employed. Actually, an example of 2-dimensional (2-D) 
case is shown in Fig. [1] (a), in which the harmonic po- 
tential is applied only along the x axis. The result of 
Fig. [1] (a) is obtained in 4-leg Hubbard model by us- 
ing the directly-extended density-matrix renormalization 
group (dex-DMRG) method [7, 8, 9] to 2-D systems. See 
Ref. [5] for other characteristic variations of 1-D system 
in different parameter ranges. 

Such specific spatial-patterns have their own novel in- 
terests, e.g., how the spin correlation structures change 
from the Mott domain to metallic edges and whether 
the holes in the metallic edges form Cooper pairs or not 
[13, [Hi. However, one clearly notices that intrinsic prop- 
erties of the Hubbard model under a fixed doping are not 
directly observable. Although the system is quite clean 
and controllable, there is no direct relationship with solid 
state physics except for a few artificial cases [12]. In 
this paper, we therefore suggest that an alternative trap, 
whose shape is box (see the right hand side panel in Fig.[T] 
(b)), enables to experimentally study the Hubbard model 
with the open boundary condition. A typical dex-DMRG 
result in the case is shown in Fig. [1] (b), which displays 
atomic density profile in the half- filled case. This result 
clearly shows that the Mott domain (^^ n^cr ~ 1) ex- 
tends over all region. By using such a trap shape, we can 
fully examine the Hubbard model under a fixed doping 
and an interaction. 

Recently, the 1-D box-shape confinement potential was 
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FIG. 1: (a) A typical dex-DMRG calculation result for the 
repulsive Hubbard model with a harmonic trap along the x 
axis as shown in the right-hand side panel. An open bound- 
ary condition is imposed perpendicular to the y axis. The 
lattice size is 20 x 4 (legs), and U/t is 10. (b) A dex-DMRG 
result for the repulsive Hubbard model with the open bound- 
ary condition for both axes. The trap shape on both axes 
is schematically depicted in the right-hand side panel. The 
filling is just the half {J2a ''T'^ = 1) and other conditions are 
the same as (a). 



actually created, and the system was successfully con- 
densed to Bose-Einstein condensate inside the box [l3]. 
According to Ref. [13], 1-D array of several small optical 
wells is also available inside the box trap in a control- 
lable manner. Thus, 1-D fermion Hubbard model with 
the open boundary condition is accessible if bosons are 
replaced by fermions. Although the box trap dimension 
is now still one, the extension to 2-D box is expected 
to be straightforward. Our expecting potential shape 
is described in Fig. [21 where 2-D box is created inside 
the x-y plane with a narrow confinement along z-axis, 
and 2-D optical lattice is loaded inside the box plane by 
adding the Gaussian walls [13] as shown in Fig. [2] or a 
standing wave light fields with shallow Gaussian envelope 
curve. Such a stage is just described by 2-D Hubbard 
model with the open boundary condition. We emphasize 
that the system becomes the best playground for simulat- 
ing controversial phenomena of High-Tc superconductors 
and other layered metal oxides. Furthermore, we would 
like to point out that ladder type models with the open 
boundary are recently good targets for advanced DMRG 
methods 0, [3 direct comparative studies are 
possible. However, it should be noted that the box shape 
is just the first step to remove the spatially dependent 
features of the filling. The temperature reduction in such 
a trap remains as the next step. 

A key issue in the repulsive 2-D Hubbard model with 
the open boundary is the stripe formation, i.e., which 
types of stripes are formed with variation of doping and 
U/t still remains controversial. Although the stripe for- 
mation have been also examined in t-J model [14], we 
focus on solely the Hubbard model in this paper. The 
reason is because t-J model requires more specific condi- 
tions in realizing its equivalent situation on atomic Fermi 
gases. Therefore, our research priority is the Hubbard 
model. In addition, we note that the Hubbard model 




FIG. 2: An expecting 2-D box-shape trap- potential with an 
optical lattice. The "endcap beams" [ij] produce Gaussian 
walls on the both ends of the x-, y-, and z-axis, while the 
optical lattice is created by additional weak endcap beams 
between the Gaussian walls. The walls along the z direction 
confines fermions much more tightly. 



requires much more efforts in numerically exploring the 
ground state than t-J model since the degree of freedom 
is much bigger than that of t-J model [15] and therefore 
the stripe formation and its profile in the ground state of 
the Hubbard model is the present-day numerically chal- 
lenging issue. 

Let us briefly review DMRG studies for the stripe for- 
mation of the Hubbard model. After an early DMRG 
work on 3-leg Hubbard ladder [iG^ , White and Scalapino 
applied the multichain algorithm of DMRG to 6-leg Hub- 
bard model and found the appearance of hole stripes 
pT^ . Afterwards, Hager et al. re-examined it in a more 
systematic manner [18]. However, their results do not 
reach to the ground state because the observed stripes 
still include the spin polarized modulation contrary to 
the Lieb-Mattis theorem as mentioned by themselves 
K [13, Although these pioneering works lacked fi- 
nal confirmations, they judged the stripe states as the 
ground state from the convergent tendencies. Thus, the 
accurate stripe profile in the ground-state are still un- 
known and U/t and filling dependence of the profile are 
also unsolved. On the other hand, the authors have re- 
cently developed dex-DMRG and confirmed that the hole 
stripe is really the ground state. The results converge 
and satisfy the Lieb-Mattis theorem. This indicates that 
the present results give the first accurate profile of the 
stripe in the Hubbard model. The extension to 2-D (lad- 
der) is achieved by parallelizing the diagonalization of 
the superblock Hamiltonian. See Ref. [7] for its method- 
ological details, performance, and accuracy. This method 
gives us a chance of direct analysis on the ground-state. 
In this paper, we confine ourselves 4-leg ladder (20 x 4) 
Hubbard-model because the present dex-DMRG easily 
reaches its ground state and enables to do systematic 
calculations varying U/t and the doping rate within our 
computational resources. When using the present dex- 
DMRG, the 4-leg case (20 x 4) almost converges within 
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FIG. 3: The dex-DMRG results of the hole distribution pro- 
files 1 — (n^,! + rii^i) for [//t = 3 (a), 5 (b), and 10 (c) in 
the Hubbard model with the open boundary condition on the 
both axes. The doped hole is just a pair (1 1, 1 1) in all cases. 
The number of sites is 20 x 4. (d) U/t dependence of the hole 
distribution along the ladder rung (y) direction at the center 
of the ladder leg (x = 10). 



3 times sweeps, which take about 3-hour s under the use 
of 128 CPU's on Altix 3700Bx2 in JAEA. We obtain the 
ground state in a doping range from p = 0.025 to 0.350 
with U/t = l^ 15. 

Let us present numerical calculation results. At first, 
we focus on a case when a hole pair (|, |) is doped into 
the half- filling. The doping rate is p{= 1 — ;^ ^ ^i^a) = 
0.025, which corresponds to a heavily underdoped region 
in High-Tc cuprate superconductors. The interest in the 
region is how the doped holes distribute, i.e., whether 
the two holes form a pair or not. If yes, then how 
the pair localizes (or delocalizes)? Such a fundamental 
question is deeply relevant to unsolved issues in heavily- 
underdoped High-Tc superconductors [19]. An alterna- 
tive experimental-stage is provided by the 2-D optical 
lattice inside a 2-D box trap. Figure [3] displays U/t de- 
pendence of the hole distribution profiles. In all cases, the 
spin-polarization completely drops to zero everywhere ac- 
cording to the Lieb-Mattis theorem. One finds for all the 
three cases (a)-(c) that the doped two holes form a stripe, 
called a hole-pair stripe. The area of the hole-pair stripe 
is much more compact than that of the single hole (| 
or I) doped case (not shown). This means that an at- 
tractive interaction effectively works between two holes. 
In addition, the dex-DMRG results show that the stripe 
shape changes from the oval to the complete stripe one 
with increasing U/t (see 2-D projected contour maps of 
(a) to (c)). This means that the shape is affected by 
the boundary edge along the ladder leg {x) direction. 
As U/t increases, the edge depression is considered to 
be defeated by the energy cost-down due to the stripe- 
formation. Fig. [3] (d), which is U/t dependence of the 
distribution profile of the hole density along the ladder 



rung {y) direction at the center {x = 10), shows that 
there is a qualitative difference between U/t = 5 and 7 
in the stripe shape. Next, in order to check the stability 
of the hole-pair stripe, we examine its leg- length depen- 
dence and confirm that the stripe keeps the same shape 
up to 40 X 4 sites model. Although the maximum limit of 
the leg- length is presently 40, any shape changes in the 
distribution profiles are not observable. Thus, one finds 
that two holes are bound and the pair sits in the most 
stable location (the center) as a stripe. 

Next, let us study interaction between the hole-pair 
stripes. The motivation comes from a question whether 
they simply merge or repel each other as multiple hole- 
pair stripes closely stay. Such a situation shows up 
when holes are sufficiently doped. If two hole-pair stripes 
merge, it means that the system prefers the hole segrega- 
tion [1^ or the stripe formation of a bi-hole pair [20] • On 
the other hand, if they are separated, then it indicates 
that the hole-pair stripe is more stable. Here, we note 
that although such a study should include the size depen- 
dence check, we concentrate on results only in a finite size 
system (20 x 4) due to lack of computational resources 
in this paper. The size dependence will be reported else- 
where [21]. Figure [4] shows U/t dependence of the hole 
distribution profile. The number of the doped holes is 8 
and the doping rate is p = 0.10, which corresponds to an 
underdoped region in High-Tc superconductors [19]. In 
Fig. m (a) (/7/t = 3), one finds 4 hole pair stripes, whose 
distance is almost equivalent. It is found in this case that 
the stripes repel each other. The surprising thing occurs 
at = 5, in which two stripes merge in the central area 
and a large stripe appears. Moreover, the merged stripe 
shows a stronger 1-D anisotropic feature (see and com- 
pare the 2-D projected contour plots), because the hole 
density does not relatively drop at the edges along the 
ladder leg {x) direction. This stripe just corresponds to a 
bi-hole pair stripe, which extends over about 4 site. Such 
a stripe structure characterized by 4 holes was predicted 
by Chang and Affleck [gO]. As U/t increases further, the 
merged stripe starts to split into two hole-pair stripes 
at U/t = 10, again, and the split becomes complete at 
U/t = 13. The bi-hole pair stripe is observable only for 
5 < U/t < 10. These behaviors indicate that the stripe 
formation occurs in not only a simple manner but also 
multiple manners depending on the doping rate and the 
interaction. The transition between these multiple struc- 
tures may be too difficult to clearly identify in solid state 
systems. Therefore, more controllable optical lattice with 
the box shape trap is promising. 

Finally, let us show a systematic result of the doping- 
rate dependence of the hole distribution profiles. The 
doping rate varies from p = 0.025 to 0.350, and U/t is 
fixed to be 10. In High-Tc superconductors, the heav- 
iest doped rate {p = 0.350) is in a overdoped range, in 
which the materials almost show metallic behaviors. Fig- 
ure [5] shows a doping rate dependence on x-directional 
profiles of hole distributions at = 2. It is found in 
Fig. [5] that there are two characteristic profiles in the 
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FIG. 4: The dex-DMRG results of the hole density profiles for 
(a) U/t = 3 (b) U/t = 5, (b) U/t = 10, and (d) U/t = 13 in 
the Hubbard model with the open boundary condition. The 
number of doped holes is 8 and the doping rate is p = 0.10. 
The other conditions are same as those in Fig. [3] 
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FIG. 5: The dex-DMRG results of the hole density profile 
along the ladder leg (x) direction of the 2nd(3rd)-leg for sev- 
eral doping. U/t is fixed to be 10. The other conditions are 
the same as those in Fig. [3] 



hole distribution. The first type seen in the relatively un- 
derdoped regime (p < 0.200(32 | 32 |)) is characterized 
by the stripe, whose modulation amplitude is relatively 
large (see Fig. [5] (a) and (b)). In the regime, the stripes 
merge and the extended ones emerge as the doped holes 
increases. The phenomenon is clearly observable from 
p = 0.075 (37T, 37i) top = 0.200 (32 j, 32 j). Inside this 



range, the bi-hole pair stripes and other stripes including 
holes more than 4 holes appear as well as the hole-pair 
ones. Especially, in p = 0.125 (35 t, 35 |), we point out 
that a bi-hole pair stripes develops quite distinctly. This 
doping rate is a characteristic one, where various singular 
behaviors stemming from the stripe formation have been 
reported in High-Tc superconductors [l^. In addition, 
one finds that two bi-hole pair stripes merge in p = 0.200 
(321, 32 1) and a wider flat plateau emerges compared to 
p = 0.150 (34 1, 34 1). The second typical profile appears 
above p = 0.250 (30 j, 30 |). The spatial modulation 
structure in the hole profile is independent on the doping 
rate. The modulations in this metallic regime are Friedel 
like oscillations [22|] rather than the stripe structures de- 
veloped due to the strong correlation [20]. 

Before closing, let us discuss finite temperature effects 
on all the present results. In atomic gases, the temper- 
ature reduction is always one of difficulties in explor- 
ing the ground-state properties of many-body interact- 
ing systems. The present results are static modulations 
of holes in the ground-state, which are expected to fluc- 
tuate in non-zero temperature. However, in sufficiently 
low temperature, the resultant tiny fluctuation does not 
matter for the static local-density probe. On the other 
hand, in relatively high temperature, the static modula- 
tions become fluctuating stripes, which can be observed 
by not local but density correlation probes [23j]. Ac- 
tually, dynamical stripe like fluctuations have been re- 
ported in High-Tc superconductors [19], in which charge 
and spin fluctuations are detectable much above the su- 
perconducting transition temperature by various probes. 
Thus, the present stripe structures can be also identified 
by using suitable correlation probes [23] in finite temper- 
ature. 

In conclusion, we suggested in atomic Fermi gases that 
the box shape trap is essential for studies on intrinsic 
properties of the Hubbard model with a fixed doping rate 
and an interaction. We applied the dex-DMRG method 
to the 4-leg Hubbard model with the open boundary 
condition in order to predict rich varieties of hole pro- 
files including stripe formations due to strong correla- 
tion. Consequently, we observed the hole-pair stripe, the 
bi-hole pair stripe, and other stripe structures includ- 
ing more than 4 holes in the underdoped region. These 
structures merge and split with the variation of U/t and 
the doping rate. On the other hand, such characteristic 
features disappear and only a monotonic oscillation pat- 
tern emerges in the overdoped one. These results can be 
directly confirmed in atomic Fermi gases, and moreover, 
more systematic studies are possible in experiments of 
atomic Fermi gases. We believe that the optical lattice 
inside the box shape trap can solve significant controver- 
sial issues like the stripe formation in High-Tc supercon- 
ductors and other metal oxides. 
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